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ABSTRACT 

This paper summarizes the results of 13 years of timing observations of a unique 
binary pulsar, PSR B1259— 63, which has a massive B2e star companion. The data 
span encompasses four complete orbits and includes the periastron passages in 1990, 
1994, 1997 and 2000. Changes in dispersion measure occurring around the 1994, 1997 
and 2000 periastrons are measured and accounted for in the timing analysis. There is 
good evidence for a small glitch in the pulsar period in 1997 August, not long after 
the 1997 periastron, and a significant frequency second derivative indicating timing 
noise. We find that spin-orbit coupling with secular changes in periastron longitude 
and projected semi-major axis (x) cannot account for the observed period variations 
over the whole data set. While fitting the data fairly well, changes in pulsar period 
parameters at each periastron seem ruled out both by X-ray observations and by the 
large apparent changes in pulsar frequency derivative. Essentially all of the systematic 
period variations are accounted for by a model consisting of the 1997 August glitch 
and step changes in x at each periastron. These changes must be due to changes in 
the orbit inclination, but we can find no plausible mechanism to account for them. It 
is possible that timing noise may mask the actual changes in orbital parameters at 
each periastron, but the good fit to the data of the x step-change model suggests that 
short-term timing noise is not significant. 
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1 INTRODUCTION 

PSR B1259— 63 was discovered in a large-scale high fre- 
quency survey of the Galactic plane (Johnston et al. 1992a). 
It is unique because it is the only known radio pulsar in orbit 
about a massive, main-sequence, B2e star (Johnston et al. 
1992b). PSR B1259-63 has a short spin period of ~48 ms 
and moderate period derivative of 2.28 x 10~^^, implying a 
characteristic age of only 330 kyr from spin-down by mag- 
netic dipole radiation. It has an orbital period of ^1237 
days, and an eccentricity of 0.87, the longest orbital period 
and largest eccentricity of any of the known binary radio 
pulsars. The inclination of the orbit to the plane of the sky 
is i ~ 36°. 

The companion, SS 2883, is a 10th magnitude star with 
a mass of about 10 and a radius of 6 Rq. Typical of 
Be stars as a class, it has a hot, tenuous polar wind and a 
cooler, high density, equatorial disk. The mass loss rate of 
the B2e star is ~ 10"® M©. Johnston et al. (1994) observed 
Ha emission lines from the disk at 20 stellar radii (R.), just 
inside the pulsar orbit of 24 R« at periastron. The density 
of the disk material is high near the stellar surface (10* — 



10^" cm~^), and falls off as a power-law with distance from 
the star. The disk is likely to be highly tilted with respect 
to the pulsar orbital plane, and PSR B1259— 63 is eclipsed 
for about 40 days as it goes behind the disk. 

Assuming the B2e star is rotating at ~70 per cent of 
its break up velocity (Porter 2001), SS 2883 then has an 
equatorial velocity ~280 kms~^. The spin- induced oblate- 
ness of the star implies an additional 1/r^ gravitational po- 
tential term in the interaction with the pulsar, known as 
quadrupole gravitational moment in typical binary system. 
This effect introduces an apsidal motion and precession of 
the orbital plane if the spin of the companion is not aligned 
with the orbit angular momentum, characterised by uo and x 
respectively, where u is the longitude of periastron passage 
and X is the projected pulsar semi-major axis (Lai, Bildsten 
& Kaspi 1995). 

In an eccentric binary system, the passage of the pulsar 
through periastron excites tidal motions on the companion 
star which then interact with the orbital motion. In the case 
of PSR J0045-7319 (Kaspi et al. 1996), tidal effects may ac- 
count for the observed evolution in binary parameters (e.g. 
Lai 1997). The tidal interactions may be enhanced by large 



© 0000 RAS 



factors if there is a resonance with oscillation modes of the 
star (Witte & Savonije 1999). For PSR J0045-7319, at pe- 
riastron the pulsar is only ~ 4 R* from the companion com- 
pared to 24 R, for PSR B1259-63. Tidal effects are strongly 
dependent on the separation of the two stars, and so may 
be insignificant in PSR B1259-63. 

Interactions between the pulsar and the disk of the Be 
star may also be important. If there is accretion from the 
disk on to the neutron star, the pulsar may be spun up or 
slowed down. For accretion to occur, the bow shock, where 
there is pressure balance between the pulsar wind and the 
Be-star disk gas, must lie inside the accretion radius, which 
depends on relative velocity of the pulsar and the accret- 
ing gas. For PSR B1259— 63, accretion is unlikely (Manch- 
ester et al. 1995; Kaspi et al. 1995; Tavani & Arons 1997), 
with X-ray observations (Hirayama et al. 1999) most eas- 
ily interpreted as originating from a bow shock well outside 
the pulsar magnetosphcrc. Disk interactions may also af- 
fect pulsar orbit directly through frictional drag. Estimates 
(Manchester et al. 1995; Wex et al. 1998) of variations in 
the Kcplerian parameters Pb, uJ and e suggest the effect is 
negligible. 

Timing models for binary pulsars have evolved over the 
years, mainly to take into account general relativistic effects 
(Blandford & Teukolsky 1976; Damour & Deruelle 1986). 
Wex (1998) extended these models to deal with long-period 
systems with main sequence companions. His model, the so- 
called MSS model, accounts for short-term periodic orbital 
effects and for the long-term secular effects caused by clas- 
sical spin-orbit coupling. 

The latest paper on timing of PSR B1259-63 (Wex 
et al. 1998) analysed observations from 1990 January to 
1996 December which covered the 1990 and 1994 perias- 
tron passages. Wex et al. (1998) were unable to derive a 
unique timing model for the system, partly because of the 
large timing noise and partly due to the long eclipses of the 
pulsed emission. They put forward two possible solutions for 
the timing. In one, timing noise dominated, its effects were 
removed through higher order derivatives of the spin fre- 
quency, and there was no orbital precession or tidal effects. 
In the second model, which they preferred, significant values 
of u! and x were obtained. These were interpreted as being 
due to orbital precession, and a large tilt between the rotar 
tion axis of the Be star and the orbital angular momentum 
vector was inferred. Two different fits were obtained for oj 
and X depending on the rmmbcr of (integer) rotations added 
during the 1994 eclipse. The preferred solution then relied 
on the physical interpretation of the values of u) and x. 

Since then, we have obtained a further six years of tim- 
ing data on this pulsar to give a total data span of 13 years, 
including two further periastron passages in 1997 and 2000. 
Our main findings are: (1) that there is strong evidence for 
a pulsar period glitch in 1997 August, about 94 days after 
the 1997 periastron, (2) that the Wex et al. (1998) timing 
solution which included orbital precession does not remain 
valid through the 1997 and 2000 pcriastrons, (3) that a solu- 
tion with jumps in i/ and i/ at each periastron, while fitting 
the data reasonably well, is unlikely to be physically viable, 
and (4) that the data are best cxijlainod by a combination 
of the 1997 August pulsar glitch and steps in projected orbit 
semi-major axis (a;) at each periastron. However, the possi- 
ble presence of timing noise and the long-duration eclipses 



conspire to make timing of the pulsar difficult. In Section 2 
we discuss details of the observations and data analysis tech- 
niques. Section 3 describes the pulse profile evolution from 
0.66 to 13.6 GHz and Section 4 describes the dispersion mea- 
sure variations observed around the 1994, 1997 and 2000 
periastrons. The timing results are described in Section 5 
with fits of the various possible models, and the possible 
interpretations are discussed in Section 6. 



2 OBSERVATIONS AND ANALYSIS 

Wc report here on observations of PSR B1259— 63 using the 
Parkes radio telescope over a time interval of nearly 5000 
days between 1990 January and 2003 June (MJD 47909 to 
52804). A total of 1031 times of arrival (TOAs), spanning 
four orbital periods, were obtained for the pulsar. The ma- 
jority of the observations were made at frequencies around 
1.4 GHz with additional observations at 0.66, 2.4, 4.8, 8.4 
and 13.6 GHz. 

At 1.4 GHz, the H-OH receiver (Thomas, Greene & 
James 1990) with a system equivalent flux density of ~40 Jy 
was used prior to 1997. From 1997, the multibeam receiver 
(Staveley-Smith et al. 1996), with a system equivalent flux 
density of ~28 Jy, was employed. Short observations with 
duration ~10 min were made at typically two-week inter- 
vals far from periastron. Longer observations with duration 
~1 hr were made on a daily basis at times close to periastron. 
The higher frequency observations were generally conducted 
close to the epochs of periastron in order to measure the dis- 
persion measure changes and the scatter broadening of the 
pulse profile (Johnston et al. 2001). The receivers at 4.8 and 
8.4 GHz had system equivalent flux densities of ~100 Jy. 
Lower frequency observations were more rarely performed. 
The combination of increased Galactic background emission 
at low frequencies, and the flat spectral index of the pulsar, 
made it hard to detect. 

The backend systems included various combinations of 
fllterbanks and correlators. Table 1 lists the details of the ob- 
serving systems. Filterbank systems were used throughout 
the observation period. At frequencies above 1 GHz, they 
generally consisted of 64 frequency channels each of width 
5 MHz (1991-1997), 96 channels of width 3 MHz (1997- 
203) or 512 channels of width 0.5 MHz (2001 onwards). The 
fllterbank backends were only capable of measuring total 
power. They one-bit sampled the data typically at intervals 
of 600 /xs for this pulsax. The data stream was written to 
tape for subsequent off-line analysis. This analysis involved 
de-dispersing and folding the data at the topocentric pe- 
riod. A digital correlator, the Caltech Fast Pulsar Timing 
Machine (FPTM; Navarro 1994), was in use between 1994 
and 2001. It had 128 MHz of bandwidth with 128 frequency 
channels. When operating in total power mode, two inde- 
pendent 128 MHz bands were available. When full Stokes 
parameters were recorded, however, only one band could 
be used. The correlator performed on-line pulse-folding and 
dedispersion to produce a pulse profile with up to 1024 time 
bins per period. Further information for earlier observations 
can also be found in Wex et al. (1998) and Johnston et al. 
(2001). 

The pulse profile of PSR B1259— 63 consists of two, 
almost equal-strength components. Each component has a 
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Table 1. Detail of the observing systems used for PSR B1259— 63. 
FB denotes a filterbank system and FPTM denotes the Caltech 

Fast Pulsar Tiiiiiiig Machiuo. 



660 MHz 



Fr6cj 


Backend. 


Year 


-D VV 


iNO. or 


Kin nf 

L\ o. or 










Channels 




0.66 


FB 


1990-1996 


32 


128 


14 




FPTM 


1995-1997 


32 


128 


13 


1.4 


FB 


1990-1997 


320 


64 


243 




FB 


1997-2003 


288 


96 


101 




FB 


2001 


256 


512 


3 




FPTM 


1994-2001 


128 


128 


444 




FPTM 


1997-1998 


64 


128 


4 


2.4 


FB 


1992-1993 


320 


64 


11 




FPTM 


2000-2001 


128 


128 


4 


4.8 


FB 


1993-1997 


320 


64 


33 




F13 


1997 


576 


192 


24 




FPTM 


1996-2001 


128 


128 


78 


8.4 


FB 


1993-1996 


320 


64 


21 




FB 


1997 


576 


192 


5 




FPTM 


1996-1997 


128 


128 


17 


13.6 


FPTM 


1995 


128 


128 


15 



Table 2. Details of the observation data spans. t+ is the first 
detection of the pulsar after the preceding periastron and r— is 
the final detection prior to the following periastron. 



Data 


Preceding 


T+ 


T — 


No. of 


span 


Periastron (MJD) 


(d) 


(d) 


TOAs 


1990.1-1990.7 






107 


18 


1990.7-1994.0 


48124 


171 


20 


187 


1994.0-1997.5 


49361 


24 


18 


443 


1997.5-2000.9 


50597 


16 


52 


237 


2000.9-2003.5 


51834 


19 




146 



steep edge and a more slowly falling edge. Manchester & 
Johnston (1995) suggested that the steep edges were the 
outside of a cone with large opening angle. In order to ob- 
tain a good value for the dispersion measure, it is important 
that the profile be properly aligned over the entire frequency 
range between 0.66 and 13.6 GHz. We chose the point mid- 
way between the mid-points of the steep outer edges as the 
fiducial point. At each frequency, a large number of pro- 
files wore summed together to form a high signal-to-noise 
template. These templates axe shown in Fig.l. To obtain an 
accurate TOA we then cross-correlated individual observa- 
tions with the relevant template. 

Table 2 lists the epoch of periastron for each of the four 
orbits so far observed. It lists the total number of TOAs 
obtained during the orbit, and indicates the first detection 
of the pulsar after perieistron and the final detection of the 
pulsar before it enters the eclipse. 

The eclipses of the pulsar, which last for ~40 days 
around periastron make the timing difficult. The problem 
is exacerbated by the fact that DM changes are observed 
over a few tens of days leading up to and immediately fol- 
lowing the eclipse due to the changing line of sight inter- 
secting the Be star wind. These DM variations need to be 
removed to obtain a good timing solution, and this involves 
making observations at several difTerent frequencies as near 




Figure 1. Standard profiles of PSR B1259— 63, summed from 
individual observations. The bar on the left side of the profiles is 
the effective resolution, including the effects of DM smearing. 

to simultaneously as possible. Multi-frequency observations 
were therefore made before and after the 1994, 1997 and 
2000 periastron passages and have been described in de- 
tail in Johnston et al. (1996) and Johnston et al. (2001). 
Generally, observations were made at 5 different frequen- 
cies, 1.2, 1.4, 1.5, 4.8 and 8.4 GHz. The three frequencies 
near 1.4 GHz were obtained simultaneously and the obser- 
vation would then be followed by one at either 4.8 or 8.4 GHz 
before again observing at the lower frequencies. 

In order to measure an accurate DM we assume that 
the DM contribution from the 1.5 kpc along the line of sight 
through the interstellar medium is 146.8 pc cm~^ (Connors 
et al. 2002). Then the extra DM contribution from the disk 
of the Be star could be computed from the relative difference 
in the timing residuals between the different frequencies. 
Typical errors using this method are ~ 0.1 pc cm^^. 

The timing properties were analysed using the pulsar 
timing program TEMPO*, which provides least-squares fit- 
ting to the pulsar rotation and orbital parameters and the 

* See http://www.atnf.csiro.au/research/pulsar/timing/tempo/ 
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dispersion measure. The MSS binary model of Wex (1998) 
was used for all fits except those in Section 5.6. To fit for 
steps in the orbital parameters, a new binary model (BTJ) 
based on the BT model (Blandford & Teukolsky 1976) was 
implemented. This allows cumulative steps in longitude of 
periastron {oj), projected semi-major axis (x), eccentricity 
(e) and binary period (Pb) to be inserted at specified times 
and also allowed setting or solving for jumps in pulsar phase 
at the specified times. 

We note that our results differ slightly from those of 
Wex et al. (1998) when analysing the same data set. This is 
because new standard profile templates were used and the 
dispersion analysis re-done for this work. 



3 PULSE PROFILE EVOLUTION 

Fig. 1 shows that the pulse profile evolves substantially over 

the observed range from 660 to 13600 MHz. At first glance 
it appears that the profile gets narrower (when plotted as in 
Fig. 1) at lower frequencies, the opposite trend to that ob- 
served in most conal profiles. However, a closer examination 
suggests that this apparent width variation is largely due 
to differing spectral index of components making up each 
peak of the profile. The outer components are stronger at 
higher frequencies whereas components on the inner edge of 
each peak become strong at lower frequencies and dominate 
the profile at frequencies below 2 GHz. The flatter spectral 
index for outer conal components is commonly observed in 
other conal pulsars (Rankin 1983; Lyne & Manchester 1988) 
and reinforces the interpretation that the steep edges define 
the boundary of a wide emission cone from a single polar 
region (cf. Manchester 1996). 



4 DM VARIATIONS 

As PSR B1259— 63 approaches periastron and moves into 
the Be-star disk, electron density variations along the line 
of sight result in DM changes as a function of orbital phase. 
Fig. 2a shows the observed variations from 1992, covering 
the 1994, 1997 and 2000 periastrons, and Fig. 2b shows ex- 
panded plots for ±75 days about the periastrons. DM vari- 
ations for the 1994 periastron were reported by Johnston et 
al. (1996), and for the 1997 periastron these were monitored 
extensively and reported by Johnston et al. (2001), while 
the 2000 variations have not been previously reported. 

At the 1994 periastron, a peak of 10.7±0.2 pc cm~^ for 
ADM was observed at r — 30; after this it decreased with 
some fluctuations until the pulsar was eclipsed. In 1997, the 
peaJi was 7.7 pc cm~^ at ~ r— 28, slightly later than in 1994, 
and again the ADM drops with some fluctuations until the 
pulsar is eclipsed at ~ r — 18. There are no multi-frequency 
data prior to the 2000 periastron. The points of ADM in 
Fig. 2 are derived by extrapolating a timing solution based 
on observations in the 60 days before these pre-periastron 
points shown in Fig. 2b. The latest pre-periastron point is 
52 days before the periastron. Based on the 1994 and 1997 
behaviour, the peak of ADM was not observed. 

The ADM fluctuations are well observed after the 1997 
and 2000 periastrons. The ADM values are much smaller 
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Figure 2. (a) Dispersion measure variations from 13 years of 
timing observation. The four vertical lines in this and subse- 
quent plots represent 1990, 1994, 1997 and 2000 periastrons for 
PSR B1259-63. (b) Expanded plots for DM variations ±75 days 
around periastron. The offset is measured relative to the value of 
146.8 pc cm-3. 



than before periastron and approach zero about 30 days af- 
ter periastron. In fact. Fig. 2 (a) shows that the ADM values 
are slightly negative in mid-orbit, suggesting the reference 
DM value of 146.8 pc cm^'' is too large. The filterbank ob- 
servations from T + 43 after the 2000 periastron to the end 
of the data set gives a mean offset of —0.2 ±0.1 pc cm~^ 
from the reference value; corresponding to an interstellar 
dispersion of 146.6±0.1 pc cm~^. 
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Figure 3. Simulated timing residuals after including jumps in v 
and V (top two panels), constant ld and x (next two panels) and 
jumps in orbital parameters (bottom four panels). Jumps were 
inserted at the second and third periastrons. See text for details. 

5 TIMING RESULTS 
5.1 Simulations 

We first used simulations to determine tlie sliape and ampli- 
tude of timing residuals corresponding to various parameters 
used in the fitting procedure. Fig. 3 shows us the timing sig- 
natures of the various parameters and the size of the resul- 
tant phase perturbations. In the figure wo have attempted 
to ensure that the resultant residuals were of the order of 
tens of milliseconds, similar to the residuals seen in the real 
data. Apart from the Co and x simulations, we introduced 
steps at the second and third periastron passages, with the 
change at the third periastron being twice the magnitude 
and of opposite sign to the change at the second perias- 
tron. Since the changes are cumulative, this has the effect 
of reversing the sign of the perturbation at the third peri- 
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Figure 4. (a) Post-fit residuals after fitting for v, v, i> and the 
Keplerian orbital parameters. The rms residual is 14.6 ms. (b) 
Post-fit residuals with the 1997 August glitch (Section 5.6) in- 
cluded, giving an rms residual is 4.3 ms. 



astron. The steps in v and v at the second periastron are 
respectively 5 x 10~^ s~^ and 5 x 10~^^ s~^. For uj and x 

we use the values from the Wex et al. (1998) model 2A of 
2 X 10""' deg yr"^ and 15 x 10"^^ respectively. The orbital 
parameter steps at the second periastron are Aw = 0?03, 
A a; = 0.01 s, Ae = 4 x 10"® and APa = 30 s respectively. 



5.2 A Glitch Near MJD 50691 

After we had been through the process of fitting the data as 
described in the subsections below, it became apparent to 
us that there was strong evidence for a glitch in the pulsar 
period near MJD 50691 (1997 August 30), about 94 days 
after the 1997 periastron. The evidence for the glitch and its 
parameters are described more fully in Section 5.6 below. In 
each of the following subsections we describe the fits both 
with and without the glitch. In all cases except the final fit 
in Section 5.6, the glitch parameters are held fixed at the 
values determined in this final fit. 



5.3 The Keplerian Timing Model 

We initially attempted to fit the entire data set using the five 
Keplerian orbital parameters and the spin frequency and its 
first two time-derivatives. Wex et al. (1998) already showed 
that this did not fit the data well. The residuals for our ex- 
tended data set are shown in Fig. 4a where the rms residual 
is 14.6 ms. The post-fit residuals after including the glitch 
shown in Fig. 4b are somewhat reduced with an rms residual 
of 4.3 ms, but still have large systematic variations. Param- 
eters from this fit are given in the second column of Table 3. 
If we add four more higher-order spin frequency derivatives, 
the rms residual reduces to ~ 2.1 ms but only at the ex- 
pense of significantly more free parameters and still leaving 
systematic residuals. We note that in the timing residual 
plots (Fig. 4-7) residual error bars arc not plotted since 
the TOA errors (typically <100 jjts) are generally smaller 
than the plotted points. In general, systematic departures 
fi-om the timing model dominate the residual plots. 
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Table 3. Four timing solutions for PSR B1259-63. 





No w, A 




V, V jumps 


Orbital jumps 


R.A. (.12000) 


13''02'"47''.65(1) 


13''02™47''.65(1) 


13'"02™47''.65(1) 


13'"02™47^65(1) 


Dec. (J2000) 


-63°50'08".7(1) 


-63°50'08".7(1) 


-63°50'08".7(1) 


-63°50'08".7(1) 


DM (pc cm-3) 


146.8 


146.8 


146.8 


146.8 


V (s-i) 


20.93692453667(7) 


20.9369245383(1) 


20.93692435(4) 


20.9369245339(9) 


i> (xlO-12 s"2) 


-0.9979077(6) 


-0.9978663(7) 


-0.9987(2) 


-0.9979383(9) 


i> (xlO-24 s-3) 


-2.11(2) 


-1.96(2) 


0.5(5) 


-1.762(7) 


Period epoch (MJD) 


50357.00 


50357.00 


50357.00 


50357.00 


X (s) 


1296.315(2) 


1296.282(2) 


1296.3264(4) 


1296.272(5) 


e 


0.8698869(5) 


0.8698832(6) 


0.8698902(2) 


0.8698872(9) 


to (MJD) 


48124.3491(1) 


48124.3494(2) 


48124.34892(3) 


48124.34911(9) 


^'b (d) 


1236.72404(3) 


1236.72360(7) 


1236.724259(7) 


1236.72432(2) 


(deg) 


138.66588(7) 


138.6644(1) 


138.66624(6) 


138.6659(1) 


Glitch epoch (MJD) 


50690.7 


50690.7 


50690.7 


50690.7(7) 


Ai/g (xlO-9 s~i) 


67 


67 


67 


67(1) 


Ai/d (xlQ-'' s-i) 


22 


22 


22 


22(1) 


0) (deg yr-i) 


— 


0.00020(1) 


— 


— 


X- (xlO^i^-) 


— 


127(5) 


— 


— 


Am,o (xlO-" s-i) 


— 


— 


0(3) 


— 


Ai>9o (x 10-15 s-2) 


— 


— 


0.8(2) 


— 


Ai/94 (xlO-'J s-1) 


— 


— 


15(3) 


— 


Ai/94 (x 10-15 s-2) 


— 


— 


0.19(5) 


— 


A!/97 (xlO-'-* s-1) 


— 


— 


-10(3) 


— 


Al>97 (xlO-l^ 






—0.30(5) 




Ai/00 (xlO-9 s-1) 






-2(3) 




Ai/00 (X 10-15 s-2) 






-0.33(5) 




Axgo (ms) 








60.3(7) 


Ax94 (ms) 








-26.3(1) 


Arr97 (ms) 








2.8(,3) 


Aajoo (ms) 








4.2(4) 


No. TOAs 


1031 


1031 


1031 


1031 


Rms residual (ms) 


4.3 


4.7 


0.78 


0.46 



5.4 Spin-orbit Coupling, Co and x 

We next attempted to obtain a timing solution by adding in 
d) and X to account for the spin-orbit coupling. As a start- 
ing point we used Model 2a from Wex et al. (1998) which 
used data up to the end of 1996. Fig. 5a shows this model 
extended to the entire data set. Clearly, its predictive power 
has failed and it does not track through the 1997 perias- 
tron. The residuals deviate immediately after the 1997 pe- 
riastron and become worse again after the 2000 periastron. 
We then allowed the orbital parameters, the spin frequency 
and its first two derivatives and Cj and x to vary; the re- 
sultant fit is shown in Fig. 5b. The residual has an rms of 
13.1 ms and shows significant deviations from white noise for 
the 1997 and 2000 orbits. The fitted values of tl) and x are 
0.00029 deg yr"! and 61 x 10"!^ respectively, both having 
changed sign from the Wex et al. (1998) solution. The post- 
fit residuals after including the glitch are shown in Fig. 5c. 
Again, this reduces the rms residual to 4.7 ms, but leaves 
systematic variations still unmodelled. Parameters from this 
fit are given in the third column of Table 3, showing that 
both uj and x remain positive. Fitting for u) and x changes 
the sign of x and only reduces the rms residual to 4.0 ms. 

5.5 Jumps in v and v at Periastrons 

Our next attempt at a timing solution involved fitting for 
frequency and frequency derivative jumps at each of the 
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Figure 5. (a) Pre-fit timing residuals using the input parameters 
from model 2A of Wex et al. (1998). (b) Post-fit residuals after al- 
lowing all parameters to vary. Wc obtain di = 0.00029(3) deg yr^ 1 
and X = 61(14) X 10^ 1^, and the rms residual is 13.1 ms. (c) 
Post-fit residuals with the 1997 August glitch included. The rms 
residual is 4.7 ms. 

four periastron epochs (in a similar fashion to that done by 
Manchester et al. (1995) on a very early data set). Fig. 6a 
shows the residuals after fitting for the five orbital parame- 
ters, the rotation frequency and its first two derivatives and 
the jumps in v and at each periastron. The rms resid- 
ual is now only 1.4 ms, and is a significant improvement on 
Fig. 4a and Fig. 5b, albeit at the expense of having 8 and 
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Figure 6. (a) Residuals after fitting for u, v, i>, the Keplerian 
orbital parameters and jumps in v and z> at each periastron, giving 
an rms residual of 1.4 ms. (b) Post-fit residuals after adding a 
glitch at MJD 50691. The rms residual is 0.78 ms. 
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Figure 7. (a) Post-fit residuals for a fit of pulsar frequency and its 
first two derivatives, the Keplerian orbital parameters and jumps 
in the projected semi-major axis, x, at each periastron, giving an 
rms residual of 11 ms. (b) Post-fit residuals after the addition of 
a glitch at MJD 50691. The rms residual is now only 0.46 ms. 



6 extra free parameters respectively. Even this fit, however, 
does not remove the systematic deviations seen between the 
periastron passages. 

The magnitude of the jumps in v arc in the range 10~® 
to 10~*s~^. For the 1990 and 2000 periastrons they are neg- 
ative, while for the 1994 and 1997 periastrons they are posi- 
tive. All the jumps in z> are of negative sign, with amplitude 
~ 10~^®s~^. However, after the ghtch is included, several of 
the jumps in v and v change sign: Ai> at the 1990 and 1994 
periastrons become positive, while Av at the 1997 perias- 
tron becomes negative. Parameters of this fit are shown in 
column four of Table 3 and the post-fit residuals are shown 
in Fig. 6b. The rms residual is now only 0.78 ms, and al- 
though some systematic variations are still seen, the fit is 
much better than for any previous fit. It should be noted 
that the jumps in z> have absorbed the overall i) term (Ta- 
ble 3). 

5.6 Jumps in Orbital Parameters at Periastron 

From the results obtained so far, it seems clear that there are 
unmodelled effects still present in the data and that these 
effects vary from orbit to orbit but are not accounted for by 
secular changes in the orbital parameters. The next logical 
step is therefore to test the effect of jumps in the orbital 
parameters, w, x, e and Pt, at each periastron. 

For the data set used by Wex et al. (1998), prior to the 
1997 periastron, an excellent fit could be obtained by fitting 
for the pulsar frequency, its first two derivatives, the Kepler 
parameters and jumps in just x at the 1990 and 1994 peri- 
astrons. However, with the full data set and fitting for Ax 
at the remaining two periastrons, large systematic residuals 
were obtained. This fit is shown in Fig. 7a which has an rms 
residual of ~11 ms. Even allowing for jumps in the other 
orbital parameters at any or all of the periastrons did not 
remove the systematic residual variations. 

Examination of the post-fit residuals from many of the 
fits showed a sharp change of gradient near MJD 50691 (1997 
August 30), about 94 days after the 1997 periastron, fol- 
lowed by a roughly exponential relaxation. This feature is 
most easily seen in Fig. 6a. Remarkably, fitting for this glitch 
in addition to jumps in just the projected semi-major axis 
at each periastron produced the almost featureless residuals 



shown in Fig. 7b. The total frequency change at the time 
of the glitch, determined to be MJD 50690.7 ± 0.7, was 
Avg = (67±1) X 10"® s-i of which Avd = (22±1) x lO"'' s"^ 
decayed exponentially with an assumed timescale of 100 
days. The final rms residual is only 0.46 ms. Only a few 
points close to the periastrons are discrepant; these could 
easily be due to errors in the DM correction. Parameters for 
the fit shown in Fig. 7b are given in column five of Table 3. 



6 DISCUSSION 

6.1 Timing Noise 

PSR B1259— 63 is a relatively young pulsar with a large 
frequency derivative which is likely to have a high level of 
timing noise (Cordes & Downs 1985; Arzoumanian et al. 
1994). Arzoumanian et al. (1994) proposed a measure of the 
timing noise to be 

A(.)=log(Mt3) (1) 

where t is the time span over which the pulsar has been 
observed. For PSR B1259-63 we have t = 10* s and for the 
measured value of the frequency and its second derivative we 
obtain As = —2.1. This is very close to the value for other 
pulsars with similar frequency derivatives, showing that PSR 
B1259— 63 is typical in this respect. 

The derived braking index, n = vv/iP', for PSR 
B1259— 63 is —36.7. Again, this merely reflects the relatively 
high level of timing noise in this young pulsar (Johnston & 
Galloway 1999). The absence of systematic residual varia- 
tions after fitting for steps in projected semi-major axis at 
each periastron (Fig. 7b) suggests that the timing noise is 
very red; higher-order derivatives in spin frequency are not 
significant. 

6.2 Steps in Rotational Parameters 

Fig. 6 shows that jumps in pulsar frequency and frequency 
derivative at periastron can account for most of the ob- 
served systematic residuals. As discussed by Manchester et 
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al. (1995), the most probable cause ol frequency jumps at 
periastron is accretion of mass and angular momentum from 
the Bc-star disk. The observed Ai/ are of order 10~^ Hz, re- 
quiring accretion of about 10^" g at the Alfven radius of 
5 X 10*^ cm. While this amount of accreted mass is not un- 
reEisonable, there axe several problems with this idea. Firstly, 
the bow shock between the pulsar wind and the circumstel- 
lar wind must lie inside the accretion radius for accretion to 
occur. Kaspi et al. (1995) and Tavani & Arons (1997) ar- 
gued that this implies unreasonably high values for the disk 
density and/or outflow velocity. Secondly, the observed 
are of different magnitudes and signs at the different perias- 
trons. While there is evidence from the unpulscd continuum 
emission that wind parameters vary from one periastron to 
the next (Connors et al. 2002), it would be surprising if the 
accretion torque were able to change sign. 

A further problem with this explanation for the ob- 
served timing residuals is that frequency jumps alone do 
not account for the observed variations; jumps in frequency 
derivative are also required. In fact, the v jumps are gen- 
erally more significant than the jumps in v (Table 3). The 
implied changes in braking torque or moment of inertia are 
large and physically improbable. These changes in v could 
result from timing noise, but the fact that they are concen- 
trated at the periastrons suggests that this is unlikely. 

6.3 Changes in Orbital Parameters 

It is clear from Section 5.4 that secular changes in uj and 
x s a sin i do not account for the observed timing behaviour 
of PSR B1259— 63, although it worked well for the first two 
periastrons. This is reinforced by the fact that the fitted 
jumps in x (Section 5.6) are of different sign at the differ- 
ent periastrons. This appears to mean that the conclusions 
drawn by Wex et al. (1998) about the inclination of the Be- 
star spin axis to the orbital plane may no longer be valid. 
We note that the x measured by Wex et al. (1998) is roughly 
consistent with the step change in x at the 1994 periastron 
described in Section 5.6. 

It is never-the-less striking that jumps at each perias- 
tron in projected semi-major axis only are required to sat- 
isfactorily model the data. Limits on changes to the other 
orbital parameters are at least an order of magnitude less 
than the values used in the simulations of Fig. 3. The sys- 
tematic residual variations seen in Fig. 4b are well modeled 
by the step changes in x at each periastron and small ad- 
justments in the other orbital parameters. If these variations 
are due to timing noise, the noise must mimic the effect of 
changes in the orbital parameters, that is, be smooth be- 
tween periastrons and rapid near periastrons, which seems 
unlikely. 

If the jumps in x are interpreted as changes in semi- 
major axis o, Kepler's Third Law implies changes in Ph of 
order thousands of seconds, orders of magnitude larger than 
permitted by the timing solutions (cf. Fig. 3). The observed 
changes in x must then be interpreted as changes in inclina^ 
tion angle i of the orbit of order a few arcsec. 

The question then is: what can cause changes in inclina- 
tion angle without making significant perturbations to the 
other orbital parameters? Furthermore, the mechanism must 
be able to produce changes of different sign and magnitude 
at different periastrons. 



There is evidence that the Be-star disk is highly inclined 
to the orbital plane both from observations of the pulsed 
emission (Johnston et al. 2001) and the unpulsed emission 
(Johnston et al. 1999), so the velocity of the disk material 
would have a significant component perpendicular to the 
orbit. In principle, frictional drag could then tilt the orbit. 
It is also known that the properties of the wind interaction 
vary substantially from orbit to orbit. The problem is that 
this effect could only change the orbit inclination with one 
sign and Manchester ct al. (1995) suggest that the fractional 
change of orbital energy due to this effect is 10"^'*, too small 
to account for the observed changes. 

The other possible mechanism appears to be tidal inter- 
actions occurring at each periastron (Papaloizou & Pringle 
1982). As mentioned in Section 1, these effects can be greatly 
enhanced by resonances with modes of stellar oscillations. 
Since there is evidently a large angle between the spin axis 
and the orbit normal, again it is possible that the dominant 
effect could be on the orbit inclination. Different relative 
phases of the tides and stellar oscillations could alter the 
effect on the orbit from one periastron to the next. Cal- 
culations to date (e.g. Witte & Savonije 1999) have only 
considered aligned systems. However, once again, the pre- 
dictions are that the effects would not be significant for PSR 
B1259— 63 because of the relatively large ratio of periastron 
distance to companion-star radius. 



7 CONCLUSIONS 

We have presented an analysis of 13 years of timing data en- 
compassing four periastron passages of PSR B1259— 63. Sig- 
nificant changes in dispersion measure occurring near each 
periastron were accounted for in the data analysis. It is clear 
that secular changes in orbit parameters due to spin-orbit 
coupling does not give a good description of the data. The 
effect of timing noise is difficult to quantify because of cou- 
pling with the orbital parameters, a problem exacerbated by 
the long-duration eclipses which occur at each periastron. 
Consequently, we cannot rule out the possibility that our 
measured parameters are affected by timing noise. However, 
excepting the glitch which occurred near MJD 50691, the 
fact that within each orbit the observed residuals show vari- 
ations with characteristic timescalc of the same order as the 
orbital period and significant shorter-term variations only 
near the periastrons strongly suggests that changes in orbital 
parameters at periastron are the dominant effect. This is re- 
inforced by the excellent fit to the data by a model contain- 
ing the MJD 50691 glitch and steps in projected semi-major 
axis at each periastron. It is clear that these steps must re- 
sult from step changes in the orbit inclination. However, we 
arc unable to offer any plausible mechanism to account for 
the observed changes. 
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